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The electronic properties of a graphene monolayer in a magnetic and a strain-induced pseudo- 
magnetic field are studied in the presence of spin-orbit interactions (SOI) that are artificially en- 
hanced, e.g., by suitable adatom deposition. For the homogeneous case, we provide analytical results 
for the Landau level eigenstates for arbitrary intrinsic and Rashba SOI, including also the Zeeman 
field. The edge states in a semi-infinite geometry are studied in the absence of the Rashba term. 
For a critical value of the magnetic field, we find a quantum phase transition separating two phases 
with spin- filtered helical edge states at the Dirac point. These phases have opposite spin current di- 
rection. We also discuss strained magnetic waveguides with inhomogeneous field profiles that allow 
for chiral snake orbits. Such waveguides are practically immune to disorder-induced backscattering, 
and the SOI provides non-trivial spin texture to these modes. 

PACS numbers: 73.22.Pr, 73.23.-b, 72.80.Vp 






cd 



T3 

O 

(N 
> 

en 

i> 

o 



X 



I. INTRODUCTION 

The physics of graphene monolayers continues to at- 
tract a lot of attention and to provide a rich source of 
interesting phenomenal— By studying the effects of the 
spin-orbit interaction (SOI) in a graphene layer, where 
symmetry allows for an "intrinsic" (A) and a "Rashba" 
(A) term in the SOI, Kane and Mele^ made a remark- 
able discovery that sparked the exciting field of topo- 
logical insulators^ For A > A/2, there is a bulk gap 
with topologically protected edge states near the bound- 
ary of the sample. This is similar to the quantum 
Hall (QH) effect but happens in a time-reversal invari- 
ant system. The resulting "quantum spin Hall" (QSH) 
edge states form a one-dimensional (ID) helical liquid, 
where right- and left-movers have opposite spin polar- 
ization and spin-independent impurity backscattering is 
strongly suppressed. The QSH state has been observed 
in HgTe quantum wells, 6 but several works 7 - - — showed 
that A is probably too small to allow for the experi- 
mental verification of this novel phase of matter in pris- 
tine graphene. Consequently, other material classes have 
been employed to demonstrate that topologically insu- 
lating behavior is indeed possible^. However, graphene 
experiment o 10 ! 11 have also demonstrated that the Rashba 
coupling A can be increased significantly by depositing 
graphene on Ni surfaces. Moreover, very recent theoret- 
ical predictions 12 suggest that already moderate indium 
or thallium adatom deposition will dramatically enhance 
A by several orders of magnitude. By using suitable 
adatoms, it is expected that in the near future both SOI 
parameters A and A can be varied over a wide range in 
experimentally accessible setups. 

In view of these developments, in this paper we study 
the electronic properties of a graphene monolayer with 
artificially enhanced SOI. Besides the SOI, we include 
piecewise constant electrostatic potentials, orbital and 
Zeeman magnetic fields, and strain-induced vector poten- 



tials. The latter cause pseudo-magnetic fields but do not 
violate time reversal invariance; for a review, see Ref. Il3l . 
While the interplay of the Rashba term A with (pseudo- 
)magnetic fields in graphene has been studied in several 
theory works beforer^~— the intrinsic SOI A did not re- 
ceive much attention so far. However, the transmission 
properties of graphene's Dirac- Weyl (DW) quasiparticles 
through barriers with arbitrary SOI have been studied 
recentl y 17 ' 18 in the absence of (pseudo-) magnetic fields. 

The structure of this article is as follows. In Sec. |TT] 
we formulate the model and construct the general solu- 
tion for piecewise constant fields. On top of the orbital 
magnetic field, we allow for arbitrary SOI parameters A 
and A, Zeeman energy 6, and we also take into account 
aspects of strain-induced fields. The homogeneous case 
is addressed in Sec. IIII1 where we determine the Landau 
level states for this problem in closed and explicit form. 
In particular, the fate of the zero modes residing at the 
Dirac point (energy E = 0) will be discussed in the pres- 
ence of the SOL Our results also apply to the case of 
a strain-induced homogeneous pseudo-magnetic fields 
Next, in Sec. II VI we study edge states near the boundary 
of a semi-infinite sample for vanishing Rashba coupling, 
A = 0. For weak magnetic fields, one then expects to have 
helical (spin-filtered) QSH edge states. Interestingly, at 
the Dirac point, upon increasing the magnetic field, we 
find that a quantum phase transition takes place between 
the QSH phase and a second QSH-like phase with spin- 
filtered edge states, considered previously by Abanin et 
al.r^ where the spin current direction is reversed. This 
spin current reversal should allow for an experimental de- 
tection of this quantum phase transition, on top of the 
obvious consequences for QH quantization rules<2°-~— In 
Sec. [Vj we turn to a mesoscopic waveguide geometry, 
where a suitable inhomogeneous magnetic field (or ex- 
change field produced by lithographically deposited fer- 
romagnetic films) defines the waveguide*^— We show 
that the SOI parameters A and A give rise to interesting 



spin texture of the resulting propagating crural states in 
such a waveguide. Finally, we conclude in Sec. I VII 



II. MODEL AND GENERAL SOLUTION 
A. Model 

Unless many-body effects are of crucial importance, 
the low-energy electronic properties of a graphene mono- 
layer are well captured by two copies of a DW Hamilto- 
nian supplemented with various terms describing SOI, 
(pseudo-) magnetic fields, and electrostatic potentials. 3 
The wavefunction corresponds to a spinor comprising 
eight components, 



/ 



*(x,y) 



^AlK 
^BiK 
^ A^K' 
^>B^K' 
^AiK' 
®BIK' 



(x,y) 



b K '{y) 



(2.1) 



The Pauli matrices (Ji =x ,y,z below act in sublattice space 
corresponding to the two carbon atoms (A/B) in the ba- 
sis of the honeycomb lattice, while Pauli matrices Si act 
in physical spin (f, I) space. Finally, the valley degree of 
freedom (K, K') corresponds to the two K points 3 - and 
Pauli matrices Ti refer to that space. Specifically, we here 
consider models where the mentioned extra terms in the 
Hamiltonian are piecewise constant along the y-direction 
and homogeneous along the x-axis. Consequently, the 
momentum p x is conserved, and we have an effectively 
ID problem in terms of the four-spinors <j> K ' K (y). The 
orbital magnetic field B z = eB (with e = ± and B > 0) is 
expressed in terms of the vector potential A.(x,y), where 
we choose the gauge 



A, 



-eB(y-c ), A y 



0. 



(2.2) 



Inclusion of the constant cq is necessary when connecting 
regions with different magnetic fields in order to make A x 
continuous. Assuming that the magnetic field is perpen- 
dicular to the graphene sheet, the Zeeman field couples to 
s z and determines the coupling constant b — g s psB/2 , 
where g s s=s 2 is the Lande factor and /is denotes the Bohr 
magneton. The full Hamiltonian then reads 3 (e > 0) 



H = vf 



O x T z [p x + ~ (A x + T z A x )j + (Ty (py + -T Z A, 



A 



Vf ~ 10 6 m/s, while the SOI couplings A and A (both 
are assumed non-negative) correspond to the intrinsic 
and Rashba terms, respectively. In wrinkled graphene 
sheets the coupling A also captures curvature effects <£ A 
constant electrostatic potential, V, has been included in 
Eq. (|2.3j) . Strain- induced forces^ 3 - lead to a renormal- 
ization of V as well as to the appearance of an effective 
vector potential, 



V + ebs z + -(cr x SyT z - cr y s x ) 



Aa z 



(2.3) 



In Eq. (|2.3[) p x = hk x is the conserved momentum in 
the s-direction, while p y = —ihd y is still an operator. 
The constant c$ in Eq. (|2.2|) can be included by shifting 
p x , and we suppose that this shift has been carried out 
in the remainder of this section. The Fermi velocity is 



JKx 



^U X y 



expressed in terms of the in-plane strain tensor u,,-, see 
Ref. [3J.The constant k can be found in Refs. Il3ll35l As 
discussed by Fogler et al.^ in many cases it is sufficient 
to consider a piecewise constant strain configuration. As- 
suming that the x-axis is oriented along the zig-zag di- 
rection, strain causes only a finite but constant A x while 
Ay = 0. This can be taken into account by simply shift- 
ing p x in this region. Below we suppose that also this 
shift has already been done. Estimates for A x in terms 
of physical quantities can be found in Refs. 13.36. The 
resulting pseudo-magnetic field then consists of <5-barriers 
at the interfaces between regions of different strain. An 
alternative situation captured by our model is given by 
a constant pseudo-magnetic field, whose practical real- 
ization has been described recently^ In that case, A x 
is formally identical to A x in Eq. (|2.2p . Unless specified 
explicitly, we consider the case of constant A x below. 



B. Symmetries 

Let us briefly comment on the symmetries of this 
Hamiltonian. In position representation, the time re- 
versal transformation is effected by the antiunitary 



operator 



37 



T = T x (-is y )C 



(2.4) 



with complex conjugation operator C and implies the re- 
lation 



TH t {k x )T- 1 = H^ t {-k x ) 



(2.5) 



for H in Eq. (J2.3I) with e = sgn(_Bj). Since H is diagonal 
in valley space, Eq. (12.5[) implies that the Hamiltonian 



H near the K' point is related to H by the relation 



rK' 



H« e (-k x )=s y [H?{k x )]*s 



(2.6) 



By solving the eigenvalue problem at the K point, we 
could thus obtain the eigenstates at K' via Eq. (|2.6[) . A 
simpler way to achieve this goal is sketched at the end of 
this subsection. 

From now on we switch to dimcnsionless quantities by 
measuring all energies in units of the cyclotron energy 
fouJc, where we define uj c = vf/Ib- The magnetic length 
£ B = (hc/2eB) 1 / 2 sets the unit of length. A field of 1 



Tesla corresponds to hui c ~ 36 meV and £b ~ 18 nm. 
Measuring B in units of Tesla, we get for the Zeeman 



coupling b = {g s fi B B/2)/hu c « 1.6 x 10- S ^B[T]. 
the dimensionless coordinate 


With 


r) = y- 2ek x 


(2.7) 


and the auxiliary quantities 




H±=E-V + b±A, u±=E-V-b±A, 


(2.8) 


we find the representation 




E-H^ +1 = 


/ v- a \ 

a* v + iX 

— iX \i+ a 

\ a) ix- j 
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(2.9) 


E-H^ = 


1 /i_ -a) 
—a /i + iA 
— iX v + —a' 
\ -a v- 


\ 
/ 




Here we introduced the standard ladder operators 




.-?+« 


a f - - - d 




(2.10) 



with [a, cr] = 1. 

According to the above discussion, eigenstates at the 
K' point for e = il could be obtained from the corre- 
sponding solutions at the K point with e = =pl. Alterna- 
tively, there is a simpler way to obtain the K 1 states as 
follows. The ID Hamiltonians H K,K (for given e) can 
be written in dimensionless notation as 



H K = 



eq 
— —Ox - %o v d n + Aa z s 

— y(7 x Sy — (7yS X ) T" J^-x^i 



ebs z , 



H 



A" 



eq 

T 

x 



i<J y d v 



Acr»s, 



+ 7r( _cr xSj/ -<T y s x ) +A x a x + ebs z 



Both Hamiltonians are therefore related by the transfor- 
mation 



H K '(A x )=(TyH K {-A x )(Ty, 



(2.11) 



without the need to invert the real magnetic field since 
this is not a time reversal transformation. As a conse- 
quence, the ID eigenstates <fi K (ry) follow from the so- 
lutions at the K point by multiplying with — ia y and 
inverting the sign of A x , 



<p K (r],A x ) = -icr y (f) K (ri,-A x ). 



C. General solution 



(2.12) 



with E — H K in Eq. (12.9[) . We construct the solution to 
Eq. (J2.13I) within a spatial region where all parameters 
(magnetic fields, strain, SOI, etc.) are constant but arbi- 
trary. This general solution will be employed in later sec- 
tions, where specific geometries are considered by match- 
ing wavefunctions in adjacent parts. Now Eq. (|2.13[) is a 
system of four coupled linear differential equations that 
admits precisely four linearly independent solutions de- 
rived in App.lAl In order to solve Eq. (|2.13[) . it is instruc- 
tive to realize that the parabolic cylinder functionS ) 38 i 39 
D p (z), obey the recurrence relations 

aD p {rj) = P D P -M, a?D p ( V ) = D p+1 ( v ), (2.14) 

with the ladder operators a,a^ in Eq. (|2.10[) . Similar 
relations for r\ — > —77 or 77 — > ir\ are given in App. [A] For 
given energy E, the order p can only take one of the two 
values 



P 



+ V- l±y/(ji + V-l) 



'I*- 



where we define [cf. Eq. (|2.8[) ] 

H = m+M- = {E-V + b) 2 - A 2 , 
v = v+v- ={E-V -bf - A 2 . 



(2.15) 



(2.16) 



For each of the two possible values for p, we then have 
two basis states, <p p and tp p , which results in four lin- 
early independent solutions. We show in App. [A] that 
the (unnormalized) solution <fi p can be chosen as 



9e=+l,p 



v_D p {-n) 

^D p (- V ) 



Xfi- 



(2.17) 



D P+ i(-v) J 
D p+ i(v) \ 

*^D p (r,) 
while if) p is taken in the form 



9e=-l,p 



( 



fl>. 



V=+l,p 



-iD_ p (-irf) 

V-D- p -i(-i-q) 



y_ij^±A D _ p _ 2 ^ l7l) ) 



(2.18) 



V>e 



-1,P 



/ i(p + l)D_ p _ a (M7) \ 
H-D- p -i{irj) 
D- p -i(irf) 



iQf-g-l) 



Next, we analyze the spatially homogeneous case. 



We now determine the spinors </> solving the DW equa- 
tion for energy E, 



(E-H K )<b(? 1 )=0, 



(2.13) 



III. HOMOGENEOUS CASE 

In this section we study an unstrained infinitely ex- 
tended graphene monolayer where the magnetic field 



B z = B (we assume e = +1) and the SOI parameters A 
and A are constant everywhere. (The electrostatic poten- 
tial V just shifts all states and is set to zero here.) We are 
thus concerned with the relativistic Landau level struc- 
ture for graphene in the presence of arbitrary SOI param- 
eters, including also the Zeeman field b. This problem 
was solved for the special case A = b = by Rashba, 16 
see also Ref. Il5l and below we reproduce and general- 
ize this solution. We focus on the K point only, since 
the spectrum and the eigenstatcs at the K' point follow 
from Eqs. (|2.6p and (J2.12I) . We also allow for a constant 
pseudo-magnetic field. When only an orbital or a strain- 
induced pseudo-magnetic field is present but not both, 
each energy level below has an additional twofold valley 
degeneracy. 

In the homogeneous case, normalizability of the spinors 
4> p [Eq. (|2.17[) ] can only be satisfied if the order p is con- 
strained to integer values p = —1,0,1,2,..., while the 
ip p [Eq. (|2.18[) ] are not normalizable. Solutions for the 
homogeneous problem thus have to be constructed using 
4> p only. Expressing the energy E (we remind the reader 
that here all energy scales are measured in units of fvjj c ) 
in terms of p [Eq. (|2.15[) ]. the sought (valley-degenerate) 
Landau levels follow as the roots of the quartic equation 

[{E + b) 2 -(p + l + A 2 )] [(25 - bf - (p + A 2 )] = 
= A 2 [(25 - A) 2 - b 2 ] . (3.1) 



For b = A = A = this recovers the standard rela- 
tivistic spin-degenerate Landau levels^ E± jTl = ±y/n for 
n = 1,2,3,... (with n = p for spin up and n = p + 1 
for spin down states), plus a spin-degenerate zero mode 
25 = (for p = 0, -1). We notice from Eq. flO) that for 
6 = 0, the combination of A and A breaks particle- hole 
symmetry, while the two couplings individually keep it. 
Furthermore, zero-energy solutions are generally not pos- 
sible except for special fine-tuned parameters. Eq. (13. 1[) 
also predicts that if E is a solution for the parame- 
ter set {p, A, A, b} then ~E is a solution for the set 
{p, A, — A,— b}. The (f> p (ri) thus represent Landau level 
states in the presence of SOI and Zeeman coupling. The 
normalization constant 1/ y/JV p , entering as a prefactor 
in Eq. (I2.17[) . can be computed analytically since D p {z) 
can be expressed in terms of Hermite functions for integer 
^Forp = 1,2,3,..., we find 



with the eigenvalue 



P 



_ V2tt P \ 
P (A/x-) 2 



(Am-)V 



(3.2) 



+^i(AV +(y- pY) + {v- pY{p + 1) 



Remarkably, for p = — 1, we find the exact normalized 
state for arbitrary system parameters, 



£ p= _i = A - b. 



(3.4) 



This unique admissible eigenstate for p = — 1 is endowed 
with full spin polarization in the 1 direction. For p = 
0, the secular equation (|3.ip becomes effectively a cubic 
equation: the solution E = A + b (i.e., v- = 0) does 
not correspond to any admissible eigenstate. The three 
allowed states are described by 



-o(v) 



\ 

A/i_^_D (-??) 

iH-vD Q (—rj) 

ivD x {-rj) 



V7^ 

V2^ 2 (l + M 2 _) + AV_*/ 2 ] 



(3.5) 



This includes a "zero-mode" partner of the p — —1 state, 
plus a pair of states obtained by mixing the spin- up n = 
and spin-down n = ±1 Landau orbitals via the Rashba 
SOI. 



A. Rashba SOI only 

For A = b = but allowing for a finite Rashba SOI 
parameter A, Eq. p. II) admits a simple solution, previ- 
ously given in Ref. |l6| and briefly summarized here for 
completeness. For p = — 1 we have the solution Q3.3p . 
which now is a zero mode, while for p = 0, 1, 2, . . ., the 
eigenenergies are given by 



Ep, a ,P — Ci 



1 + A 2 



(3.6) 




1/2 



with a, (3 = ±. According to our discussion above, here 
Eq ,± t — = should be counted only once, with eigen- 
state^ _ ex (0,D (—r]),0,—i\Di(—r])), while ■Eo,±,+ = 
±Vl + A 2 correspond to a particle/hole pair of first Lan- 
dau levels modified by the Rashba SOI, with eigenstates 
<t>o,±,+ ex (0,XD (-7 1 ),±iVTTX I D a (-r 1 ),iD 1 (-7 1 )). We 
thus get precisely two zero-energy states. 
For small A, we find the expansion 



E P -i,±,+ = ±(1 
E Pt ± t _ = ±(1 



-A 2 /2)^ + 0(A 
-A 2 /2)^ + 0(A 



which shows that the states E Pl ± t + and E p+ ±,±^, which 
form a degenerate Landau level for A = 0, are split by a 
finite A. 



0-1 fa) = 



(2tt) 1 /4 









D (-v) 



(3.3) 



B. Intrinsic SOI only 

Let us next consider the case A = 0, where one has a 
QSH phase 4 for B — and A^0. Now the Hamiltonian 



is block diagonal in spin space and the eigenstates be- 
come quite simple even for finite Zeeman coupling, since 

we can effectively work with the bi-spinors <fi^ \ (y) for 
spin s =t / 1= ±. We easily obtain the (unnormalized) 
eigenstates with p E No in the formic 



k ( n _ ( fp,±, 8 D p ^i(-r]) 
-Dp(-V) 

Vp,±,sD p -l(-T)) 



(3.7) 



%±,M 



where the eigenenergies follow from Eq. ()3.1[) . 



E Pt ±, s = sb±y/p + A 2 



(3.8) 



We employ the notation 

u p>±>a = E Pt ± >a - E 0> _ s>s = ±Vp + A 2 - sA. (3.9) 

For p — 0, the second index in P ,±, S and -B P ,±, S should 
be replaced by — s, i.e., there is only one solution for given 
spin (and valley). Note that -Eo,+4 hi the present nota- 
tion corresponds 40 to the solution (|3.3|) . When 6 = 0, 
interestingly enough, A does not lift the spin degeneracy 
of the Landau levels except for the zero mode (p = 0)4i 
A Zeeman term with 6 = A restores a true doubly- 
degenerate zero-energy state for p = again. In Sec. IIVI 
we show that this implies a quantum phase transition. 



ing expressions are not illuminating and too lengthy to 
be quoted here. Only the p = — 1 state in Eq. (|3.3[) re- 
mains exact for arbitrary parameters. We here specify 
the leading perturbative corrections around the special 
cases above, and then show the generic behavior in two 
figures. 



Expanding around the Rashba limit of Sec. IIII Al 
which is justified for b, A <C 1, we get the lowest- 
order perturbative correction to the finite-energy (i.e., 
p y^ 0, — 1) Landau levels Q3.6[) in the form 



$E Pt ± t+ = —8E Pt ±- 



(A 2 A + 6) 



V(l + A 2 ) 2 +4pA 5 



(3.10) 



Expanding instead around the intrinsic SOI limit of 
Sec. IIII Bl we find the following small- A corrections to the 
Landau levels in Eq. p.8p : 40 For p — 0, the state E 0<+ ^ 
corresponding to the exact solution (|3.3[) is not changed 
by A to any order, while -Eo.-.-f obtains the lowest-order 
correction 



SE , 



2(A - 6)A 2 
46(6- A) + 1' 



C. General case 

Although the quartic equation (|3. 1 [) can be solved an- 
alytically when both SOI couplings are finite, the result- 



J 



The corresponding eigenstate is, however, not a spin- 
t state anymore. For p > 0, the eigenenergy E p> ± >s 
[Eq. (|3.8[) ] acquires the perturbative correction 



8E V 



±,a 



±- 



s\ 2 



2(A-s6)(AtvV+A 5 ) 



2^/^+A 2 



l + 46fs6± v/pTA 2 
I 



(3.11) 



We now consider two different SOI parameter sets 
consistent with the estimates in Ref. [El and show the 
complete evolution of the Landau levels from the weak- 
to the strong-field limit. In Fig. [1] numerical results 
for the few lowest-energy Landau levels are depicted for 
A > A/2, corresponding to a QSH phase for B = 0. 
The (valley-degenerate) spin-split levels corresponding to 
the A = A = 6 = zero mode exhibit a zero-energy 
crossing at B m 11 T for the chosen SOI parameters. 
This crossing signals a quantum phase transition from 
the QSH phase, which survives for sufficiently small B 
and A > A/2, to a peculiar QH phase for large B. As we 
discuss in Sec. IIVI one then again has helical edge states 20 
but with reversed spin current. Similar crossings can oc- 



cur for higher Landau states as well, as is shown in Fig. [5] 
for a parameter set with A < A/2 where no QSH physics 
is expected. For even larger B not displayed in Fig. [21 
we find an E = crossing where the Rashba-dominated 
small-S phase turns into the helical QH phase. 



D. Spin polarization 



Given the Landau level eigenstates, it is straightfor- 
ward to compute the spin-polarization densities Si(y) = 



^ 



-1- 



A=0.65meV 
1=0.15 meV 



10 
B[T] 



15 



20 



FIG. 1: (Color online) Low-lying Landau level energies (in 
units of the cyclotron energy hw c ) vs magnetic field B (in 
Tesla) for the SOI parameters A = 0.65 meV and A = 
0.15 meV. For small B, this corresponds to the QSH phase, 
A > A/2. For better visibility, the deviation from the respec- 
tive A = A = b = level has been magnified by a factor 10 
for each curve. 



tq 



A=1.5 meV 
X=6.5 meV 




-1 



10 
B[T] 



15 



20 



FIG. 2: (Color online) Same as in Fig.[Qbut for A = 1.5 meV 
and A = 6.5 meV. 



#+f # (i = x, y, z). We find S x (y) = 0, while 



S y (y) 



v — p 



pD p ^iDp 



-D P D P+1 , (3.12 



2AC, 



V 2 Dl_ x 



/<- 



S-^*L)Dl 






■D 



P+i 



where D p = D p {—rj). In the absence of the Rashba term 
(A = 0), the in-plane component S y vanishes identically, 
since then the eigenstates are simultaneously eigenstates 
of s z . For finite A, integration over y yields a vanishing 



expectation value for the overall in-plane polarization, 
but the Rashba coupling still induces local in-plane spin 
polarization. The case A = b = has been discussed in 
detail by Rashba. 16 



IV. QH EDGE STATES FOR INTRINSIC SOI 

In this section, we consider the edge states correspond- 
ing to the relativistic Landau level problem in Sec. IIIII 
when a boundary at y = is present. We focus on 
the case of purely intrinsic SOI, A = 0, but the physics 
should be qualitatively unchanged for A < A. In the 
region y < we then have a homogeneous magnetic field 
B z = +B, i.e, e = +1. (For a pseudo-magnetic field, this 
holds at the K point while at the K 1 point, B z — > —B z .) 

Since the problem of edge states in graphene has been 
studied extensively before, some remarks are in order at 
this point. In fact, putting A = b = A = 0, our results 
are consistent with those of Refs. [22iJ23||42h44| reporting 
chiral QH edge states in graphene. On the other hand, 
the B = model is equivalent to the continuum limit 
of the Kane-Mele model 4 and thus exhibits helical QSH 
edge stateSf 6 - (The helical state has a pair of counterprop- 
agating ID modes with opposite spin polarization.) The 
Kane-Mele model with (A, b) ^ but without orbital 
magnetic field has recently been studied^ and a quan- 
tum phase transition from a (generalized) QSH phase for 
b < A to a quantum anomalous Hall (QAH) phase for 
b > A has been predicted. It is worthwhile to stress 
that the QSH effect survives even when time-reversal 
symmetry is broken. In the QAH phase, one has chi- 
ral edge states moving in the same direction for both 
spin polarizations^ 6 - The valley analogue of this quan- 
tum phase transition has also been studied4£ Further- 
more, for the 2D topological insulator realized in HgTe 
quantum well structures, a related transition has been 
predicted— by including the orbital field but omitting 
the Zeeman term. 

However, the Zeeman term is crucial in graphene near 
the Dirac point: for A = and b ^ 0, spin-filtered helical 
edge states (similar to the QSH case) emerge again.] 20 i 49 
Our results below show that this QSH-like phase is sep- 
arated from the "true" QSH phase by a quantum phase 
transition at b = A. Albeit both phases have spin- filtered 
edge states, they differ in the direction of the spin cur- 
rent. This feature should allow to experimentally dis- 
tinguish both phases and to identify the quantum phase 
transition separating them. In practice, one may reach 
this transition simply by changing the magnetic field. 

Normalizability of the wavefunctions for y — > — oo 
implies^, that the only allowed solutions follow from 
the 4> p spinors in Eq. (|2.17|) . while the tp p solutions 
[Eq. (|2.18p ] have to be discarded. Since we do not 
have to impose normalizability at y — > oo, the order p 
is not constrained to integer values and can now take 
any real value consistent with suitable boundary condi- 
tions at y = 0. For given conserved momentum k x and 



spin s, the solutions for p yield the edge state spectrum, 
E s (k x ). Note that for finite magnetic field and k x < 0, 
the distance from the boundary is set by \k x \. Putting 
A = 0, possible solutions <f> '± s (y) must be of the form in 
Eq. (|377J}, with energ y E P! ±J given by Eq. $£%})■ While 
p € No in Sec. IIIIB1 we now consider arbitrary real p. 
To make progress, we have to specify boundary condi- 
tions at y = 0. We investigate two widely used boundary 
conditions, namely the zig-zag edge and the armchair 

edgc&amsi 



A. Zig-zag edge 

For a zig-zag edge with the last row of carbon atoms 
residing on, say, sublattice A, the microscopic wavefunc- 
tion must vanish on the next row outside the sample, 
belonging to sublattice B. In the continuum limit, since 
the x-axis here points in the zig-zag direction, the lower 
component of the spinor <j>p± s [Eq. (|3.7|l ] has to vanish 



at y = 0j i For both spin directions s 
the condition 



±, this yields 



D p {2k x ) = 0, (4.1) 

which has to be solved for the energy, expressed in terms 
of p as E s — sb ± \Jp + A 2 . At the other Dirac point, 
the lower component of the spinor (f>p± s should vanish 
at y = 0, where Eq. (|2.12j) implies the condition 

^,±, s iVi(2M = 0, ( 4 - 2 ) 

with fp,±,s m Eq. (|3.9[) . It is not possible to find simul- 



taneous solutions to both Eqs. (|4.1[) and (J4.2I) . Possible 
states are thus confined to a single valley: the bound- 
ary condition does not mix the valleys but lifts the KK' 
degeneracy. Remarkably, for s = ± and arbitrary k x , 
Eq. (|4.2p is satisfied by the K' solution for p = in 
Sec. IIIIB1 with E s {k x ) — s(b — A), i.e., we find a pair 
of "flat" states. For all other states, Eq. (|4.2p simplifies 
to condition (|4.1|) with p -> p — 1 (and if — > A''). We 
mention in passing that for A = this condition reduces 
to Eq. (9) in Ref. |44J- Equation (|4.ip can be solved in 
closed form for k x — > — oo using asymptotic properties of 
the parabolic cylinder function. To exponential accuracy, 
with n G No we find 

lot- |2n+l 

p = n+ |2 7L e ~ 2fc *. (4.3) 

Numerical analysis of the above equations recovers the 
expected spin-filtered helical edge states^ - for 6 > A, but 
the continuum approach used in this paper fails to give 
clear evidence for the helical QSH edge states for b < A. 
As pointed out in Ref. 49, under the zig-zag boundary 
condition one needs a more microscopic description in 
order to capture these states. The "flat" states above 
are remnants of the sought QSH edge states, but the 
continuum model is not sufficient to describe their proper 
dispersion relation. We therefore turn to the armchair 
boundary condition. 



B. Armchair edge 



Under the armchair boundary condition, we instead 
impose ^ + &%' = and *§ + \&f = at the bound- 
ary, with "J in Eq. (|2.1[) . This boundary condition mixes 
the valleys and involves both sublattices. Since in our 
coordinate system the x-axis is parallel to the zig-zag 
direction, we first rotate the system by tt/2 and then im- 
pose the boundary condition at y = 0. Written in the 
original coordinates, we find (for each spin direction s) 



,± s Up—iyZfc x ) zh Up{Zfc x 



0. 



(4.4) 



We note that the relative phase between the K and K' 
components is not fixed by the Dirac equation, which 
is diagonal in valley space. However, the only relative 
phase compatible with the boundary condition imposed 
simultaneously on both sublattices is ±1. Each of the two 
conditions in Eq. (|4.4I) may thus be imposed separately. 
We have checked that the numerical solution of Eq. (|4.4j) 
for A = recovers the known results for the QH edge 
state spectrumi 20 i 44 In addition, for B = 0, the armchair 
edge is known 50 ' 52 to yield QSH edge states. 

Our numerical results for the dispersion relation 
E s ,±{k x ) for the armchair edge are shown in Fig.[3j where 
± corresponds to the symmetric or antisymmetric lin- 
ear combination in Eq. (|4.4[) and the magnetic field is 
B = 15 T. The main panel shows results for A = 6 meV. 
Then A > b, and we have the (generalized) QSH phase. 
Indeed, for E = we find the helical edge state, where 
the right- (left-)mover has spin s =t (s =1)- The in- 
set of Fig. [3] is for A = 0.3 meV, where A < b and the 
spin-filtered helical QH phased is found. Here we have 
spin s =4- ( s =t) f° r the right- (left-)mover. Hence the 
spin current differs in sign for A > 6 and A < 6, with 
a quantum phase transition at A = b separating both 
phases. This feature should allow for an experimentally 
observable signature of the transition. 



V. SPIN STRUCTURE IN MAGNETIC 
WAVEGUIDES 

In this section, a spatially inhomogeneous situation is 
considered, where a magnetic waveguide 2 ^— along the 
x-direction can be realized. Since the problem remains 
homogeneous along the x-direction, p x = hk x is still con- 
served. For the physics described below, the Zeeman 
coupling b gives only tiny corrections 2 ^ and will be ne- 
glected. Moreover, there are no valley-mixing terms such 
that we can focus on a single valley. 

We distinguish a central strip of width 2L (the "waveg- 
uide"), —L<y<L, and two outer regions y < —L and 
y > L. In the central strip, we shall allow for arbitrary 
SOI parameters A and A. In addition, strain may cause 
a constant contribution to the vector potential, A x , and 
a scalar potential, V. The magnetic field in the cen- 
tral strip is denoted by B c . For \y\ > L, we assume 



In- • * « i 



3 



» 8 8 



8 ! 





. 1 ' ' ' » 1 _ 






0.02 


— A A - 




A * A 





.....' A 


« » 4 * * | • 




o • 


0.02 


_ o • _ 








, • 1 - 



a © 6 6 



• « * ♦ ♦ i . I ° 



i s= - , antisym. 

o s= - . sym. 

* s=+ , antisym. 

• s=+ , 5>'m. 



jt Z n 



FIG. 3: (Color online) Dispersion relation E s ,±(k x ) of a semi- 
infinite graphene sheet with an armchair edge at y = 0, ob- 
tained numerically from Eq. (14. 4[) . We use A = 0, B = 15 T, 
A = 6 meV, and the + ( — ) sign is for the symmetric (anti- 
symmetric) valley combination in Eq. (|4.4[1 . Inset: Same for 
A = 0.3 meV. 



that all strain- or SOI-related effects can be neglected, 
A = A = A x = V = 0. In principle, by lithographic 
deposition of adatoms, one may realize this configura- 
tion experimentally. For y < —L, the magnetic field is 
B z — B > 0, while for y > L, we set B z = eB, where 
e = 1 (e = — 1) corresponds to the parallel (antiparal- 
lel) field orientation on both sides. For e = — 1, we take 
B c = 0, while for e = +1, we set B c = —B. 

The setup with e = — 1 could be realized by using a 
"folded" geometry , 53 ! 54 cf. recent experimental studies^ 5 - 
Note that when the magnetic field changes sign, one en- 
counters "snake orbits," which have been experimentally 
observed in graphene pn junctions. 56 For the e = — 1 con- 
figuration, we have uni-directional snake orbits mainly lo- 
calized along the waveguide, while for e = +1, we get two 
counterpropagating snake states centered near y = ±L. 
For A = A = A x — V = 0, both cases (e = ±1) have been 
studied in detail in Ref. [27|. Technically, one determines 
the eigenstates and the spectrum, E(k x ), by matching 
the wave functions in the three different regions, which 
results in an energy quantization condition. This method 
can be straightforwardly extended to the more complex 
situation studied here by employing the general solution 
in Sec. |IT]for the central strip. 

Before turning to results, we briefly summarize the pa- 
rameter values chosen in numerical calculations. We take 
a magnetic field value B — 0.2 T, and the waveguide 
width is 2L = \/8£b ~ 40 nm. The strain-induced pa- 
rameters in the central strip are taken as A x = — 16^m _1 
and V — — 20 meV. These values have been estimated for 
a folded setup, 5 - where V comes from the deformation 
potential. We consider two different parameter choices 
for the SOI couplings: Set (A) has A = 13 meV and 
A = 3 meV, corresponding to the QSH phase. For set 
(B), we exchange both values, i.e., A = 3 meV and 



A = 13 meV. 



A. Antiparallel case: Snake orbit 

Let us first discuss the e = — 1 configuration, where the 
magnetic field B z differs in sign in the regions y < —L 
and y > L. The dispersion relation of typical low- 
energy ID waveguide modes is shown in Fig. [4] For 
k x — > — oo the centers of the quantum states are lo- 
cated deep in the left and right magnetic regions, far 
from the waveguide. Thus one has doubly-degenerate 
dispersionless "bulk" Landau states. With increasing k x 
these states are seen to split up. The dominant splitting, 
which is already present for A = A = 0, comes from the 
splitting of symmetric and anti-symmetric linear combi- 
nations of the Landau states for y < —L and y > L with 
increasing overlap in the waveguide region^ Asymptot- 
ically, the dispersion relation of all positive-energy snake 
states is E{k x —$■ +00) ~ hvpk x ^- For intermediate k x 
and (A, A) 7^ 0, however, we get spin-split snake states 
out of the previously spin-degenerate states. The spin 
splitting is mainly caused by the Rashba coupling A and 
disappears for A — > 0, cf. the inset of Fig. |4j 

The zero-energy bulk Landau state (for k x — > —00) 
shows rich and interesting behavior in this setup. While 
for k x — > +00, we expect one pair of snake states with 
positive slope and one pair with negative slope, for the 
studied parameter set and range of k x L, there is just one 
state with negative slope while three branches first move 
down and then have a positive slope. Accordingly, at the 
Dirac point (E = 0), Fig. [4] shows that there are three 
right-movers with different Fermi momenta and different 
spin texture. Two of those states are indicated by stars 
(*) in the main panel of Fig.2]and their local spin texture 
is shown in Fig. [5] Evidently, they are mainly localized 
inside the waveguide and have antiparallel spin polariza- 
tion. We find spin densities with S x — for both states. 
For the Rashba-dominated situation in Fig. [SJ spin is po- 
larized perpendicular to the current direction and has a 
rather complex spatial profile. 



B. Parallel configuration 

Next we come to the e = +1 configuration, where the 
magnetic field is +B for \y\ > L and —B for \y\ < L. One 
therefore expects two counterpropagating snake states in 
the a;-direction localized around y = ±L. The corre- 
sponding spectrum is shown in Fig. [5] We focus on pa- 
rameter set (B), since for set (A), the spin splitting is 
minimal and less interesting. The spectrum consists of 
two qualitatively different states, namely states of bulk 
Landau character for large |A; X |L, and a set of propagat- 
ing waveguide modes.— The spectral asymmetry seen in 
Fig. [S] for all propagating modes, E(—k x ) 7^ E(k x ), is 
caused by the strain (.A^) -induced shift of k x . Such a 
spectral asymmetry may give rise to interesting chiral- 




FIG. 4: Dispersion relation of the lowest few energy branches 
for a strained magnetic waveguide with e = — 1 and SOI in 
the central strip of width 2L. Energies are given in units of 
Hvf/L. The main panel is for parameter set (B). The stars 
refer to the states further studied in Fig. \S\ Inset: Same for 
set (A). (See main text for details.) 




FIG. 6: Same as Fig. [4] but for the setup with e = +1 and 

parameter set (B). Solid and dashed curves are for better 

visibility only. The two states indicated by stars are studied 
in Fig. 
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FIG. 5: (Color online) Spin density profile Sy, z (in arbitrary 
units) vs y/L for the two E = right-moving states indicated 
by stars in the main panel of Fig. [4] The left star corresponds 
to k x L = 1.745, the right star to k x L — 2.629. Inset: Particle 
density, p, and current density, J x (which is the only non- 
vanishing component), in arbitrary units vs y/L. We show 
the result only for k x L — 1.745, since k x L — 2.629 yields 
practically the same. 



ity and magnetoasymmetry effects^ The spin texture 
is shown in Fig. [7] for a pair of right- and left-moving 
states with E = 1.2hvp/L, cf. the stars in Fig. [6] We 
observe from the main panel in Fig.[7]that the spin polar- 
ization of both states is approximately antiparallel. Be- 
cause of their spatial separation and the opposite spin 



direction, elastic disorder backscattering between these 
countcrpropagating snake modes should be very strongly 
suppressed. The inset of Fig. [7] shows the current den- 
sity profile across the waveguide. Although the profile is 
quite complex, we observe that the current has opposite 
sign for both modes. 



VI. CONCLUDING REMARKS 

In this work, we have studied the magnetoelectronic 
properties of monolayer graphene in the presence of 
strong intrinsic and Rashba-type spin-orbit couplings. 
According to a recent proposal, 12 large intrinsic cou- 
plings may be realized by suitable adatom deposition 
on graphene. We have presented an exact solution for 
the Landau level states for arbitrary SOI parameters. 
When the intrinsic SOI dominates, by increasing the 
magnetic field, we predict a quantum phase transition 
from the quantum spin Hall phase to a helical quantum 
Hall phase at the Dirac point. In both phases, one has 
spin-filtered edge states but with opposite spin current 
direction. Thus the transition could be detected by mea- 
suring the spin current either in a transport experiment 
(e.g., along the lines of Ref. \5&t) or via a magneto-optical 
experiment. 

In inhomogeneous magnetic fields, especially when also 
strain- induced pseudo- magnetic fields are present, inter- 
esting waveguides can be envisioned. Such setups allow 
for snake states, where spin-orbit couplings result in a 
spin splitting. In a double-snake setup, there is a pair 
of counterpropagating snake states that carry (approx- 
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We here show only the case e = +1 near the K point; all 
other cases follow analogously. Solving the first and last 
equations for fa and fa, respectively, we find 




1 



V— 



a) fa 



FIG. 7: (Color online) Spin density S VjZ (in arbitrary units) 
vs y/L for the two states indicated by stars in Fig. [6] The left 
(right) star corresponds to a left- (right-)mover with k x L = 
— 1.46 (k x L = 2.83). Note that the spin polarizations of both 
states are approximately antiparallel. Inset: Particle current 
profile J x (in arbitrary units) vs y/L for both states. Black 
solid curve: k x L — 2.83. Dashed red curve: k x L = —1.46. 



imately) opposite spin polarization. This implies that 
scattering by elastic impurities is drastically suppressed. 
The resulting spin textures can in principle be detected 
by spin resolved ARPES (see, e.g., Refs. [lfl and [591) or 
spin-polarized STM measurements. 

We hope that our predictions can soon be tested ex- 
perimentally. 



For \X- = 0, one has the solution (|3.3[) instead, while for 
V- = there are no solutions. The second and third 
equations then yield two coupled second-order ordinary 
differential equations for fa and fa, 

[a* a — v)fa — iXv-fa = 0, 
iXfi-fa + (a T a + 1 — m)03 = 0. 

Solving for fa yields 



-(a'a - v)fa 



iXv_ 



and we thus arrive at the equation 



T><j)2 — [(a' a + 1 — n)(a/a — v) — A/i-^-] fa = 0. 



Since the operator T> commutes with the "number oper- 
ator" a'' a, the sought solutions for (j>2 span the kernel of 
T> where the (f>2 are eigenstates of a^a, 



a^a fa, p {rf) = p fa, P (r)). 
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Appendix A: Derivation of the eigenstates 

Here we provide some details concerning the derivation 
of Eq. (J2.17I) ; the notation below is explained in Sec. [TTJ 
First, additional relations like Eq. (|2.14[) can be stated, 



aD p (-ri) 
aD p {iri) 



-pDp-i{-r)), a t L» p (-r/) = -D p+1 (-rj), 
-iD p+1 (ir)), a 1 D p (ii]) = -ipDp^ir/) . 



We wish to construct the solution <fi — 
satisfying Eq. (j2~l~3"]) . 

/ 







V- 


a 
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(<h 


at 


v + iX 







h 





—iX fi + 


a 




h 





a* 


V- 


) 


\fa 



which implies the two solutions in Eq. (|2.15[) . With a^a — 
(?y/2) 2 — 1/2 — d 2 /dr] 2 , the eigenvalue equation for fa is 
just the differential equation of the parabolic cylinder 
functions. 38 



drj 2 



Mv) = o, 



which has the four (linearly dependent) solutions 
{Dp(n), D p (-?]), D-p-^irj), £>_ p _i(-i?7)}. Given the so- 
lution for fa, all other components in follow by using 
the recurrence relations of the D p functions, see, e.g., 
Eq. (|2.14p . After straightforward but lengthy algebra, 
we obtain the four solutions (also quoted for e = —1) 
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I 

V 



-i{v-p) 
Xfi— 



D P+ i(v) J 



-iD_ r 



-p\-irj) 
y_D_ p _ 1 (-irj) 
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\ / pD p _ x (-rj) \ 

v_D p (-r]) 

\ / -iD_ p (ir)) 

V-D- p -i(ir)) 



\ 



Q-p)(p+i) 

A/i 



H-D p (-t)) 



-j(^-p-l) jj 

\V- P~ 



MJ 



D- p - 2 (iv) J 



I -i(p + l)Z>_ p (-Mj) \ 



/'- 



_D 



»(m-p-i) 

A 



-p-1 



(-irj) 



D-p-i(-iv) 



d p+i(v) 
li-D p {rj) 

iAi ^Dp-M) 



H_D_ p _ x {iri) 

^^D_ p _ t (i V ) 



x 

\ A^_ 



For a given energy E, Eq. (|2.13|) admits precisely four 
linearly independent solutions for <f>. However, Eq. (I2.15[) 
implies two possible values for p, i.e., we have the freedom 
to choose just two out of the four quoted eigenstates (for 



given e) and then allow both values of p in Eq. (|2.15l) . 
Our conventions for these two basis states are specified 
in Eqs. (|237|) and (J2~T%1) in the main text. Thereby we 
have obtained all possible solutions to Eq. (|2. 13|) . 
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